The Warm Ionized Medium in the Milky Way and Other Galaxies 

John S. Mathisf] 

Received ; accepted 



x Dept. of Astronomy, Univ. of Wisconsin; 475 N. Charter St., Madison WI 53706, USA; 
mathis@uwast.astro.wisc.edu 



- 2 - 



ABSTRACT 

There are now observations of several emission lines from the "Warm Ionized 
Medium" (WIM) (or, equivalently, the "Diffuse Ionized Gas") of the local ISM, 
the Perseus Arm in the Milky Way, and also in several other galaxies. Interesting 
features of these observations include the great strength of [N II] A6563 (~ Ha in 
some cases) and the fact that [S II]A6717/[N II] A6583 is almost the same (~0.6 
- 0.7) in all locations and objects. Other line ratios (e.g., [O III]A5007/rL3) vary 
considerably. 

This paper presents simple photoionization models that reproduce the observed 
spectra, providing extra heating beyond that supplied by photoionization is 
assumed (see Reynolds, Haffner, & Tufte 1999). The same extra heating was used 
for models of all stellar temperatures being combined together, although it could 
easily depend on 

With observed gas-phase abundances (not solar), the line ratios in the local 
arm at b = 0° are fitted with no extra heating and (S/H) = 13 ppm, as opposed 
to solar (~ 20 ppm). Local gas observed at b = —35° requires extra heating of 
about T_25 = 0.75, where T_ 2 5 is the extra heating in units of 10~ 25 erg H _1 s -1 . 
In the Perseus Arm there are similar results: little extra heating is required at 
\z\ = 500 pc, and T_ 2 5 = 3.0 is needed at \z\ = 1.2 kpc. To fit the observations, 
the gas-phase composition in the Perseus Arm must be reduced as required by 
the Galactic abundance gradient observed for H II regions. The requirements for 
NGC 891 (the best observed other galaxy) at \z\ = 1 kpc and 2 kpc are similar 
to the Perseus Arm: little or no extra heating near the plane (1 kpc in this case) 
and T_25 ~ 3 at \z\ =2 kpc. In NGC 891 there is also an increase of A5007/Ho; 
with \z\ that can only come about if most of the ionizing radiation is supplied by 
very hot stars (type 04: T* ~ 50000 K). Either their radiation must propagate 
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from the plane to high \z\ through very little intervening matter, or else the stars 
are located at high \z\. The total power requirement of the extra heating is £ 
15% of the power to photoionize the WIM without extra heating. 

Extra heating enhances [O II]A3727/H/3. Figure 1 shows that there is a spread 
in the predicted values, but the ratio can serve as a useful diagnostic of extra 
heating. The [S III] A9065,9531 lines (see Figure 2)are not useful in diagnosing 
extra heating. 
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1. Introduction 

It is well known that very diffuse ionized gas is a major component of the interstellar 
medium (ISM) of the Milky Way (Reynolds 1991, 1993). In some (but not all) other 
galaxies about 40% of the Ha is diffuse (Rand 1996; Hoopes, Walterbos, & Greenawalt 1996; 
Domgorgen & Dettmar 1997; Greenawalt et al. 1998; Otte & Dettmar 1999), as opposed to 
H II regions (although it is not completely clear that the definition of WIM is the same for 
extragalactic objects as for the Galaxy.) The diffuse gas will be referred to as the "Warm 
Ionized Medium" (WIM) in this paper, although it is also called "the Reynolds layer" for 
the Milky Way gas and also the "diffuse ionized gas." The properties of the WIM, especially 
regarding its heating, are reviewed by Reynolds, Haffner, & Tufte (2000). 

The source(s) of the ionization of the WIM have been debated, but photoionization by 
radiation from OB stars seems to be the only viable source for the bulk of the material. Even 
for the local Milky Way, the power required for the ionizations is very large: 10 - 15% of the 
energy produced by OB stars, or all of the total power from supernovae. It is not obvious 
how the ionizing radiation from OB stars can penetrate the diffuse H I that is observed in 
all directions from the Sun. However, it seems likely that OB stars ionize enough of the 
neutral material surrounding them to allow the required fraction of ionizing radiation to 
escape (Miller & Cox 1993; Dove, Shull, & Ferrera 2000). Since the Sun does not possess this 
ionizing power, there are apparently no directions in our sky that allow ionizing radiation 
produced by the Sun to escape to great heights from the Galactic plane. 

The hypothesis that the WIM is photoionized by hot massive stars is supported by 
the agreement of its spectrum with models of very dilute H II regions (e.g. Domgorgen 
& Mathis 1994; Greenawalt, Walterbos, & Braun 1997; Martin 1997; Wang, Heckman, & 
Lehnert 1998; Sembach et al. 2000). The main reason for the general agreement is that the 
observed spectrum is heavily biased to low stages of ionization ([N II], [S II]), with weak lines 
of the higher stages that are abundant in ordinary H II regions ([O III], [S III]). This bias 
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is a natural feature of nebulae photoionized by comparatively dilute radiation fields. The 
models predict the observed anticorrelation of [N II] /Ha and [S II] /Ha with the intensity 
of Ha. However, [O I]A6300 is weak in the only two Galactic WIM regions for which it has 
been measured (Reynolds 1985, Reynolds et al 1998), both in the galactic plane. This line 
is strong only at the outer edges of H II regions, where there is both 0° (with abundance 
locked to H° by resonant charge exchange) and warm free electrons. Present models of the 
local WIM (e.g., Sembach et al. 2000) assume that some of the ionizing radiation escapes 
from the Galaxy, so the resulting A6300 is not detected. The model reflects this leakage by 
being matter bounded, with the fraction of neutral H, X(H°), ~0.1 at the outer edge instead 
of ~1. 

The heating of the WIM is less well understood than the ionization. Reynolds & Cox 
(1992) suggested that some form of heating in addition to that resulting from photoionization 
might be required. At present there are several new observations that make the case for extra 
heating very strong (Reynolds, Haffher, & Tufte 1999), and suggest that the extra heating, 
relative to the photoelectric heating, decreases with density. Thus, photoionization heating 
dominates in dense regions, such as near the midplane or in bright H II regions, but extra 
heating may become stronger than photoionization heating in the very diffuse WIM. 

The present paper will interpret several sets of observations in a general way. The 
emphasis will be to compare observations relatively close to the plane of the respective 
locations with those farther from the plane. Photoionization models will be constructed 
with simple geometries (i.e., volume averages of uniformly dense nebulae.) Combining 
the photoionization models from assumed individual stars in order to mimic the actual 
complicated situation will result in conclusions that are plausibility arguments rather than 
detailed models of the real situations. Firm conclusions can be reached as regarding heating 
requirements beyond photoionization. More elaborate modeling seems unjustified in view of 
the geometrical complexity, such as fractal structure, that is found in the actual ISM. 
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The new observations that I have chosen to model are discussed in §2 and listed in Table 
1. The general characteristics of the photoionization models are discussed in §3, and results 
regarding individual regions are in §4. Discussion are given in §5 and a summary in §6. 

2. Observations 

Older observations of the Milky Way WIM within the local spiral arm have been 
summarized and modeled in Sembach et al. (2000). Their models are close to those of 
Domgorgen & Mathis (1994). This paper will discuss the lines [S II]A6717, hereafter "[S II]"; 
[N II]A6583, hereafter "[N II]"; [O III]A5007, hereafter "[O III]"; [O I]A6300, hereafter 
"[O I]", and the quantity tpue = (He + /He)/(H + /H), the fractional ionization of He relative 
to H, derived from He I A5876/H/3. I will fit line ratios [N II]/Ha, [N II]/ [S II], [O I]/Ho;, 
[O III] /Ha, and the abundance ratio -?/>He in various objects. This paper addresses the 
[S II] /[N II] ratio, while previous papers (Domgorgen & Mathis 1994; Sembach et al. 2000) 
considered the absolute values of the lines. 

Excepting the local WIM, where there is good spatial resolution, ratios from near galactic 
midplanes are not suitable for modeling because of contamination by relatively dense H II 
regions, including radiation scattered by dust. The emission becomes too faint for reliable 
measurements at large distances from the plane. 

New observations of the local WIM and for the Perseus Arm (Haffner, Reynolds, & Tufte 
1999) were taken with the Wisconsin Ha Mapper (WHAM) near the Galactic plane and also 
down to b = —35°. The Perseus Arm WIM can be identified by the velocity shift produced 
by galactic rotation, and I will consider only \z\ = 500 pc and 1200 pc in the Perseus Arm 
in order to avoid the relatively dense H II regions near b = 0°. Haffner et al. (1999) found 
[N II] /[S II] -0.6 in the local ISM at all b. 

Reynolds et al. (1998) have measured [O I] in the local Galactic plane. In the Perseus 
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Arm they observed it at two radial velocities with strengths [O I] /Ha = 0.044, and also 
<0.012, at \z\ = 330 pc. I adopted [O I] /Ha = 0.028 at \z\ = 500 pc as a reasonable 
compromise. There are as yet no measurements of [O III] for local gas at b = —35° and for 
the Perseus Arm. 

While there are many observations of various galaxies, I choose to concentrate on 
NGC 891 because it has the greatest variety of observed lines of any galaxy. However, we 
must recognize that the observations used here are from very deep spectra (Rand 1997, 1998) 
taken at one slit position perpendicular to the plane at the location of an especially bright 
filament. The physical conditions at the location at a given \z\ were not typical of the entire 
galaxy. The average of the line ratios at the same \z\ on the opposite sides of the galaxy was 
used. The asymmetry above and below the plane is not large, so the average is physically 
meaningful. 

The adopted observations and other quantities explained below are given in Table 1 
for each of the regions that are modeled. The adopted observed value is an average of a 
substantial number of individual values near a particular value of z. 

There are common properties of the spectra of the local WIM, the Perseus Arm, and 
the WIM in NGC 891 that must be explained. One of the most striking is the constancy of 
[S II] /[N II] with \z\. It is readily explained by a constant (S/N) plus extra heating (Reynolds 
et al. 1999). However, without extra heating a strong variation in the (S/H) with \z\, highly 
correlated with a different strong correlation of (N/H), is needed (Rand 1998; Reynolds et al. 
1999). The reason is that the [N II] and [S II] are not produced in precisely the same region: 
S + is much more confined to the outer regions of the Stromgren sphere than is N + . 

There are substantial differences among the three regions. The local WIM and, by 
assumption, the Perseus Arm have modest [O I] and ^He (Reynolds et al. 1998; Reynolds 
& Tufte 1995), while NGC 891 has appreciable [O I], larger ^He at high \z\, and, especially, 
large measured [O III]/H/5. 
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3. Modeling the Observations 

My approach was to compute individual photoionization model H II regions, each 
with a single central star and a specific value for the "ionization parameter," U, which 
uniquely determines the ionization structure of the model. The parameter U is defined in 
terms of a uniform density, n e , by U = L/(AnR 2 s cn e ), where Rs is the Stromgren radius 
and L is the H-ionizing photon luminosity of the exciting star. I use a related quantity, 
q = L 50 f 2 riH, where / is the filling factor and L 50 = 10~ 50 L(photons s _1 ). It is easy to show 
that U = K q 1/3 , with K = [10 50 o;|/(367r(l + y)c 3 )] 1/3 , where y is the He/H ratio by number 
and the recombination coefficient of H to levels n > 2. In order to vary q, I varied / and 
assumed a constant L and n(H). The ionizing spectra were based on the atmosphere models 
of Kurucz (1991). The nebular models were calculated with the code used in Domgorgen & 
Mathis, with atomic constants updated to those in Pradhan & Peng (1995). 

It is important to realize that each model is characterized by the mean ionization 
parameter, q, but contains within it a whole range of variations of the local ratio of 
(photon/electron) densities, since the stellar radiation field becomes weaker near outer 
portions of the model. One cannot associate a particular measure of %,hoton/^e, such as might 
be inferred from observations of the emission measure and distance from ionizing sources, 
with the model that has that same average value of n p hoton/?V However, volume-averaged 
models probably mimic the complicated physical situation of the real ISM as well as any 
assumed radial variation of density, or even a guess regarding the three dimensional structure 
of the gas density. The actual ISM is clumped on scales covering several orders of magnitude, 
and the radiative transfer through it is very poorly represented by any current models. 
A volume average seems to be a reasonable average. The assumption of uniform density 
becomes serious for integrations along particular lines of sight piercing the model. 

Most models were assumed to have a spatially uniform mass density and rate of extra 
heating. This case represents a constant deposition of heating per H nucleus, independent of 
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the state of ionization or other physical conditions. Each volume-averaged model shows how 
uniform extra heating per H nucleus affects the averaged spectrum. The results cannot be 
used to determine the dependence of the extra heating rate on n e , as has been considered 
in the Perseus Arm by Reynolds et al. (1999). Such a prediction would involve detailed 
assumptions regarding the radial geometry of the nebula around the exciting star(s). 

The final fit to the observations was determined by combining two or three individual 
H II region models with possibly different values of q. The agreement of the combination 
with the observations was made by minimizing the statistic 



where i refers to the quantity to be fitted (such as [N II]/Ho;,) k to individual H II region 
models in the combination, and the are discussed below. The Obsi are the observed ratios; 
Pk,i is the volume-averaged ratio in each model. The minimization was accomplished by 
adjusting the a^, the fraction of the ionizing photons contributed by the individual model k. 
Only three individual models were considered for any combination, but usually only two were 
satisfactory. In view of the arbitrary nature of the geometry of the individual H II region 
models (see below), including more individual models in the combination would not have 
increased physical understanding. 

In order to assess the quality of fit of models, I assign "errors," to each quantity 
to be fitted. This "error" is not an assessment of the actual true uncertainty, but rather 
a fitting parameter that reflects how closely I require the final superposition of individual 
photoionization models to conform to the adopted observational values. The e» are listed in 
Table 1 for each ratio in each position. 

Another approach to fitting the observations might be to estimate the numbers of 
stars of various spectral types within clusters that are exciting the WIM, to combine their 
luminosities, and then to compute individual photoionized models with that stellar input. 
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This procedure is somewhat arbitrary as well. Sometimes the stars in a cluster act as a 
collective source of radiation, as for a supershell (e.g., Dove et al. 2000), but often stars 
produce their own individual H II regions that contribute to the WIM. The ionizing radiation 
from very young clusters is likely to be absorbed by the local gas, producing a rather bright 
H II region, and contributing to the WIM through leakage of some of its ionizing radiation. 
When the observed intensity is along a path spanning an entire galaxy, the actual geometry 
is very complex. 

Using rays piercing the outer regions of models might be just as appropriate as volume 
averages, especially at large \z\ if the exciting stars are found at smaller \z\. This case 
represents looking through the outer edges of the ionized region. I investigated using this 
type of averaging, and it sometimes achieved excellent fits to the observations. However, the 
results presented here are all volume averages because the geometry of real nebulae seems 
so convoluted that volume averages seem more appropriate. A major problem with using 
integrations along paths is that the radial density distribution within the model must be 
assumed arbitrarily. Mathis and Rosa (1991) showed that volume averages of models with 
varying q can well represent the spectra of various regions within individual H II regions. 

Other input parameters are the abundances of He, C, N, O, Ne, and S, all relative to H. 
Values for each region are discussed in §[| below. The abundance of Ne has a minimal effect 
on the observed ratios: a 10% increase in Ne decreases [O III] by 1.7% and the other ratios 
by < 1.2%. Carbon has almost no effect on the results because it is not a major coolant in 
the ionized regions, as opposed to the neutral ISM, which lacks the powerful coolants [O II], 
[N II], and recombination cooling. If (S/H) is increased, with the other abundances held 
constant, [N II] decreases slightly because of the extra cooling. [S II] and [S III] increase 
as (S/H) - 8 , resulting in [S II]/ [N II] oc (S/H). For a similar reason, [S II]/ [N II] is almost 
inversely proportional to (N/H). 

The fraction of neutral H at the outer edge of the model, X e dge(H°), must also be 
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assumed. Radiation-bounded models have X edge (H°) = 0.9, at which point the nebular 
temperature, if there is no extra heating, is decreasing rapidly. Continuing to X edg0 (H°) = 
0.95 gives almost identical results, but is probably no closer to reality than X e d g e(H ) = 0.9 
because of the neglect of the dynamical effects that are present near ionization fronts. I also 
considered gas-bounded models with X edge (H°) = 0.1, representing leakage of some ionizing 
radiation to larger \z\. The fraction of escaping ionizing photons depends on q, since low 
values of q imply a very gradual transition at the outer edge of the model. For the range 
of q I considered (—2.6 > log(g) > —4.8), the fraction of escaping photons is ~ 20 - 55%, 
respectively, when X e dgc(H°) = 0.1. 

The rate of extra heating per H nucleus, in units of 10~ 25 erg H" 1 s _1 , is described by 
the quantity T_25- The extra heating is assumed to take place only out to X(H°) = 0.5. The 
same value of T_ 2 5 was used for models of all stellar temperatures being combined together, 
although T_25 could easily depend on T*. The heating and power requirements are discussed 
under the various regions whose spectra were fitted (§§). The power required is considered in 
§5. 

The parameters that have major effects on the results at a given value of q are T*, 
X e d g c(H°), r_ 25 , and, marginally, (O/H). The effects follow intuition. Increasing T* shifts 
the ionization balance to higher ions. Even low values of q show measurable [O III] (~ 0.1 
H/3) if T* = 50,000 K, while no models with T* = 35,000 K produce appreciable [O III]. 
Decreasing X e d g e(H°) strongly decreases the emissions from low stages of ionization, especially 
[O I]A6300, which is strongly concentrated to the outer edges of the models. For instance, 
models with X e d gG (H°) =0.1 have A6300 a factor of ~5 weaker than those with X edge (H ) 
= 0.9. Increasing (O/H) lowers the general level of forbidden lines by providing extra 
cooling through [O II]A3727. Extra heating increases all emissions by increasing the nebular 
temperature. However, [O I]A6300/[S II] A6717 is almost independent of extra heating because 
the line emissivities depend upon temperature-dependent collision strengths and ionization 
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fractions, despite the somewhat different Boltzmann factors. 

Dust is a source of uncertainty that does not affect the results significantly. WIM models 
have much lower values of q than bright H II regions, implying a larger fraction of neutral 
H at a particular fraction of R s . Thus, in WIM models the opacity contributed by H is 
relatively large. At the ionization edge of H, the absorption cross section of an H atom is 
6.3 lCr 18 cm 2 and of dust is ~ 3 lCT 21 cm~ 2 IT 1 (e.g., Mathis, Rumpl, & Nordsieck 1977). 
Thus, if the fraction of H° is > 5 10~ 3 , as it is over most of the volume for these low values of 
q, the opacity of H° is more than 10 times the dust opacity. 

4. Results 

In this section we discuss the constraints placed on the physical nature of the successful 
combinations for local gas, the Perseus Arm, and for NGC 891. 

4.1. The Local ISM 

For CNO abundances in the local ISM I ignored solar values and used the observed 
gas-phase (O/H) = 319 per 10 6 H (Meyer, Jura, & Cardelli 1998); (N/H) = 75 ppm (Meyer, 
Cardelli, & Sofia 1997); (C/H) = 140 ppm (Sofia, Fitzpatrick, & Meyer 1998). These give a 
larger (N/O) than solar, presumably because an appreciable amount of O, but little or no N, 
is incorporated in silicates and oxides in grains. For the local ISM, the assumed abundances 
are very similar to those called "Warm Neutral Medium" in Sembach et al. (2000). The 
(S/H) required to fit the WIM spectrum at b = 0° is 13 ppm, so (S/N) = 0.17. Solar 
(S/N) = 0.22 ± 0.08 (Grevesse, Noels, & Sauval 1996). The local abundances at b = -35° 
were assumed to be the same as in the plane because it seems very unlikely that significant 
amounts of N are incorporated into dust grains (Whittet et al. 1998; Tielens at al. 1999). 
Unless N and S can be released by destruction of grains above the galactic plane, it seems 
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unlikely that their gas phase abundances will increase with \z\. 

For 6 = 0°, several combinations fitted without extra heating. The best fit, matching 
the observations virtually perfectly was produced by 56% of the ionizing photons from 
models with T, = 35,000 K, (spectral type 09 V; Vacca, Garmany, & Shull 1996), 12% 
from 40,000 K (spectral type 07 V) with X cdgc (H°) = 0.1, and the rest from T* = 45,000 K 
(type 05.5). The predictions of the combination are listed in Table 1. The relatively high 
stellar temperature required probably represents the hardening of the stellar radiation that 
leaks from denser H II regions because of the absorption of the softer photons. Other good 
combinations used as little as 32% of the ionizing photons from T* = 35,000 K, with the 
other ionizing photons divided between T* = 40,000 K radiation bounded and gas bounded 
models. The only problem with these models is their prediction that A5007/H/3 ~ 0.08, while 
the observed is 0.18. Models with = 45,000 K fitted all of the adopted observations almost 
perfectly. For any particular combination there is a range of possible values of q. 

Even a small amount of extra heating at b = 0° is problematic with the abundances 
adopted here. The major challenge for the models mentioned above is to avoid overproducing 
[N II] in the plane, and extra heating makes the problem worse. Heating as low as T_ 2 5 = 
0.4 increases [N II] /Ha to 0.65, while the observations give 0.60 ± 0.02 (Haffher et al. 1999). 
Reducing the assumed (N/H) is not an option because the ratio is observed (Meyer et al. 
1997). 

For the local arm (-y(LSR) ~0 km s" 1 ), the spectrum at b = —35° differs from that at b 
= 0° by having [N II]/Ha ~ 0.80 in place of 0.60 at b = 0° and [S II]/Ha increased by about 
the same ratio. There are no measurements of [O I] or [O III]. Extra heating is required to 
produce the stronger [N II] /Ha and [S II] /Ha. A combination of models similar to those for 
b = 0°, but with T_25 = 0.75, fits the observations: 50% of the ionizing photons from T* = 
35,000 K, and the remainder from 45,000 K. The predictions of the combination are listed in 
Table 1. If we assume that the [O I], [O III], and ipu c are the same at b = —35° as for 0°, 
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we can eliminate models with T_ 2 5 = 0. They produce either too little [N II] or else much 
too much He I, with the best fit being an uncomfortable accommodation of the two errors, 
predicting ?/>He ~0.63. The upper limit in the plane is ipn e < 0.4. Extra heating at the level 
of T_25 = 0.7 can produce an excellent fit to all ratios, including those assumed to be similar 
to the plane ([O I], ipu e , and [O III]). Observations of ^He, [O I], and [O III] would be very 
useful. 

We have seen that a combination of T* = 35,000 K and hotter stars produces good fits 
to the spectrum of the local gas. Domgorgen & Mathis (1994) used a single T* = 38,000 K; 
Sembach et al. (2000) used 35,000 K. 

4.2. The Perseus Arm 

Abundances in the Perseus Arm were estimated by fitting the observations of the WIM 
rather than using the galactic abundance gradient. They were (C, N, O, Ne, S)/H = (85, 45, 
250, 71, 8) ppm, respectively. If they had been assumed to be the local gas phase abundances, 
reduced by a factor proportional to Dq 066 (Dq is the distance from the Galactic center), they 
would have been (85, 50, 213, 85, and 9.3) ppm, respect ivelyQ. In view of the assumptions 
within the models, the differences between the two sets of abundances are not significant. 

Like the results for the local spiral arm, several combinations of models fit the 
observations in the Perseus Arm at \z\ = 500 pc. No extra heating is required. Examples 
of successful combinations are 55% of ionizing photons from T* = 35,000 K and the rest 
from 40,000 K. The predictions (see Table 1) are excellent. Another very good fit arises 

2 The exponent 0.066 is the mean of those found in Afflerbach, Churchwell, & Werner (1997) 
for O, N, and S, very similar to the gradient in B star abundances found by Gummersbach et 
al. (1998). The distance to the Perseus Arm from the Sun was assumed to be 2.5 kpc and the 
observations to the Perseus Arm were at I = 145°. The reduction factor in this case is 0.72. 
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from 35,000 K (70%) combined with 45,000 K, A edge (H°) = 0.1. The best fit is with 66% 
of the ionizing photons from models with T* = 35,000 K, 4% from 50,000 K, X edg0 (H°) = 
0.9, and 30% from 50,000 K, X c d ge (H°) = 0.1. However, such very hot stars do not seem 
appropriate for the Perseus Arm. The range of reasonably successful combinations illustrates 
the difficulty of being very specific about the ionizing spectrum. 

At \z\ = 1200 pc, there are no observations of either [O I] or [O III], but [N II]/Ha ~1 
and [S II] /Ha = 0.6. Models without extra heating fail to produce enough [N II]; even a T* 
= 50,000 K model produces [N II]/Ha of only 0.75. Increasing (N/O) and (S/O) to provide 
the required [N II] and [S II] is not attractive because the abundances would have to greater 
than solar, while the (O/H) would need to kept low in order to keep the nebula hot. 

Models with T* = 35,000 K (50% of photons) combined with T* = 45,000 K fits the 
observations at \z\ = 1200 pc well (see Table 1) with extra heating of r_25 = 3. The models 
produce relatively large amounts of [N II] and [S II], the opposite of those at 500 pc. Models 
with lower densities in the outer regions also produce an excellent fit. 

The power requirements of the extra heating will be discussed in §|. 

4.3. NGC 891 

The major difference between NGC 891 and the Galaxy is that in NGC 891, A5007/H/? 
increases with \z\ to ~0.75 at \z\ = 2 kpc, as opposed to <0.2 in the plane of the Milky Way. 
Other differences are the relative large strength of [O I] /Ha (~0.1 in NGC 891 vs. < 0.03) 
and ipHe (0.55 in NGC 891, vs. 0.4 in the Perseus Arm). 

In order to avoid the emission from dense H II regions near the plane, the smallest 
\z\ I considered in NGC 891 is 1 kpc. The Perseus Arm at that \z\ requires extra heating. 
However, NGC 891 has about the same [N II] /Ha at 1 kpc as in the plane of the Milky Way, 
so I assumed no extra heating in NGC 891 at \z\ — 1 kpc. Even low values of T-25 (~ 0.4) 
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make it difficult to obtain enough [0 I] without overproducing [N II] and [S II] . 

The abundances in NGC 891 were assumed to be (C, N, 0, Ne, S)/H = (85, 45, 250, 71, 
8) ppm, respectively. The Perseus Arm is the same except it used (O/H) = 213 ppm. The N, 
O, and S were found by fitting the spectra at both \z\ — 1 and 2 kpc, with no extra heating 
at \z\ — 1 kpc. The C and Ne have little influence on the spectra and were assumed to be 
the same relative to N as in the Galaxy. 

It is not possible for local gas phase abundances, (N, 0)/H = (75, 320) ppm, to fit the 
rather low [N II] /Ha in NGC 891 at the same time as the ipn e and [O III]. With this high an 
(N/H), [N II] is overproduced by the hot stars needed for [O III]. Assuming appreciable extra 
heating in NGC 891 at \z\ = 1 kpc makes all of these problems worse. Such extra heating 
cannot be ruled out if the compositions are taken to be free parameters. We would require 
NGC 891 to have (N/H) < 45 ppm (about half solar) when integrated across the galaxy, and 
a gas phase (N/O) < 70% of the local value. 

Good fits for NGC 891, \z\ = 1 kpc are found with ~50% of the photons from T* = 
35,000 K and the rest from T* = 45,000 or 50,000 K (Types 05.5 or earlier). The cool stars 
are needed to reduce the ?/>hc to ~0.5, and the hot stars in order to produce enough [N II] 
and ^He- Using T* = 40,000 K (type 07.5) for the hot stars does not fit because its [O III] 
is too weak. If stars with T* > 45,000 K are included, the major problem in NGC 891 is 
not in fitting the [N II], as it was for the Perseus Arm, but in producing enough [O I]. The 
best combination is a composite of three models: T* = 35,000 K (48%), T* = 50,000 K, 
A edge (H ) = 0.9 (32%), and T* = 50,000 K, X cdgc (H°) = 0.1 (20%). It produced an excellent 
fit (x 2 /^ — 0.2, where v is the number of line ratios to be fitted = 5.) Table 1 lists its 
predictions. However, the range of combinations with (say) x 2 / u < 0.6 is large. There are 
successful combinations with values of q ranging over an order of magnitude, all with the T* 
> 45,000 K for the hot stars. 

At \z\ = 2 kpc for NGC 891 we assume the same abundances as for \z\ = 1 kpc, which 
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is an important constraint because the spectra are rather different: [N II] and [S II] ~ Ha 
instead of 0.6 Ha at 1 kpc, [O III] = 0.75 H/3 as opposed to 0.55, and tpu e = 0.85 vs 0.55. 

Because of ip~Re, ~90% of the ionizing radiation must arise from very hot stars (T* 
~45, 000-50, 000 K) instead of ~ 50% at 1 kpc. To produce high ^ He , higher values of T* are 
required in the WIM than in H II regions because the H and He have much large fractions 
of neutrals at any particular fraction of the nebular radius. These ions compete for + 
ionizing photons (hv > 35.5 eV) much more successfully than in the more compact nebulae. 
In addition, strong extra heating is required at 2 kpc in NGC 891. Extra heating rate T_ 2 5 
= 3 is the least that provided a good fit to the spectrum; it required 14% from T* = 35,000 
K and the remainder from T* = 50,000 K (73% radiation bounded; 13% from gas bounded.) 
The hot stellar temperature could be reduced if more extra heating (T_ 2 5 = 4) is allowed. 
With r_ 25 = 4.5, T* = 45,000 K can be used in place of 50,000 K. 

The presence of very hot stars at large \z\ in NGC 891 is not unreasonable. There is a 
great deal of dust found well out of the plane (Howk & Savage 1997), and in other edge-on 
spirals as well (Howk & Savage 1999). If there is dust, there can be star formation and 
associated hot stars found nearby. 

Must these very hot stars be at large \z\, if the WIM at large \z\ is photoionized by 
starlight? In order to simulate "chimneys" with hot stars near the plane, I computed models 
with central voids of O.6-R5. At a given [N II]/Ha, the volume-averaged [O IIIJ/H/3 is about 
the same or only slightly weaker than in the case of uniform density. However, these models 
depend on having no gas between the exciting stars and the outer regions. For low values 
of q, H° and He efficiently absorb photons with hv > 35.5 eV, the ionization edge of +2 . 
Unless there is almost a void in the gas below 2 kpc and exciting stars near the plane, there 
would be so little +2 that extra heating will increase [N II] and [S II], but [O III] would be 
too weak. 
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4.4. Predictions of [O II] and [S III] Lines 

The strongest line from the WIM is almost surely [0 II]A3727 because of the observed 
weakness of [0 III], the abundance of oxygen, and the efficiency of [0 II] as a coolant. Figure 
1 summarizes the predictions of [O IIJA3727/H/3 of the models of various stellar temperatures. 
The solid line is for % = 50,000 K; the short dashes are for 45,000 K; the dotted is for 40,000 
K; the long dashes are for 37,500 K; the long/short dashes are for 35,000 K. At various points 
along each line, error bars indicate the range of line strengths among the various values of q 
considered for each model. The lowest value of A3727/H/9 corresponds to the lowest log(g) 
(= —4.8), but the highest occurs at various values around log(g) ~ —3.5. 

The composition of the models shown by lines and error bars were (C, N, O, Ne, and S) 
= (85, 45, 250, 71, and 8) ppm, respectively, which is the composition used for NGC 891. 

The points in Figure 1 show the mean A3727/H/3 for T* = 50,000 K (circles), 40,000 
K (squares), and 35,000 K (triangles), for nebulae containing more heavy elements than 
those used for the lines in the figure: (C, N, O, Ne, and S) = (141, 75, 320,118, 13) ppm, 
respectively. These are the abundances adopted in this paper for the local WIM. The 
difference in the points and lines shows that A3727 is influenced by the cooling from other 
ions, while the rough agreement between the two compositions shows that A3727 is a major 
coolant, especially if there are large amounts of extra heating. Our composite models for the 
WIM used roughly equal amounts of T* = 35,000 K and 50,000 K, which is not far from T* 
= 40,000 K as regards A3727/H/3. 

From Figure 1 we see that there is a general correlation of A3727 with heating, but 
with a large scatter, especially at low extra heating rates. The A3727/H/3 cannot be used to 
estimate low heating rates accurately, but it can distinguish between T_ 2 5 = and r_ 25 = 
3. From observations of A3727/H/5, Figure 1 can be used to estimate the mean amount of 
extra heating in the WIM. The figure cannot be used to determine the density dependence 
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of the extra heating, since the assumed geometry of the models does not correspond to the 
unknown density distribution and clumpy structure of the real WIM. 

Figure 2 is similar to Figure 1, except that it shows [S III]A(9065+9531)/Ha for the same 
models as in Figure 1. The values are almost proportional to (S/H) for large rates of extra 
heating. We see that [S III] is not a good diagnostic of physical conditions within the WIM. 

5. Discussion 

In spite of significant differences among the spectra of the WIM in the local ISM, the 
Perseus Arm, and NGC 891, there a considerable similarity in the simple photoionization 
models fitting them. No extra heating is required, and little tolerated, at low \z\, where 
[N II] /Ha = 0.5 - 0.6. At larger \z\, where [N II]/Ha ~1, T_ 2 5 ~ 3 is needed in at least the 
outer regions of the plasma where N and S are largely singly ionized. There is little evidence 
from these photoionization models regarding the relative amounts of extra heating in the 
regions near the exciting star(s), where the photoionization heating dominates the extra 
heating if it is assumed to be uniform. 

At large \z\ in the Perseus Arm and NGC 891, very hot stars (T* > 45,000 K) are needed 
to explain the comparatively large [N II] /Ha, [S II]/Ha, and [O I] /Ha ratios observed. These 
high stellar temperatures might represent the hardening of the radiation escaping from denser 
H II regions closer to the plane, since the softer photons are absorbed first by the H and He. 
In NGC 891very hot stars, and not hardening of the stellar radiation, are needed to explain 
the [O III] observed at \z\ = 2 kpc. For the low values of q required for the WIM, the gas 
producing the hardening would absorb the + ionizing photons required for the A5007 line. 
These hot stars could be in the plane of the galaxy, powering a superbubble (Dove et al. 
2000) if there is very little gas between them and \z\ — 2 kpc. 

These combinations of models fail to predict precisely what values of the ionization 
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parameter, or q, are required to produce the observed spectra. For any particular object there 
is a rather wide range of values that are satisfactory. The other geometrical uncertainties, 
such as using volume integrals, simplistic stellar atmospheres, or requiring the same heating 
for all stellar temperatures at a given location, make predictions of the dilution of the 
radiation field quite speculative. 

The propagation of the ionizing photons from the hot stars to the WIM may be the most 
serious objection to the idea the photoionization is the dominant process operating in the 
WIM, since only a very small amount of neutral H is completely opaque. The fractal nature 
of the ISM is shown in molecular clouds (e.g., Elmegreen & Falgarone 1996; Stutzki et al. 
1998) and in H I (Green 1993). The percolation of ionizing radiation through hierarchical 
density structures has not yet been investigated quantitatively at an adequate level, but 
probably ionizing photons from an O star can penetrate large distances through the lowest 
density regions within the fractal ISM. 

While photoionization from stars seems to be very plausible for providing most of the 
ionization of the WIM, it is also very likely that other sources contribute. Candidates 
are shocks (Raymond 1992) caused by supernovae, as well as turbulent mixing layers 
(Slavin, Shull, & Begelman 1993) produced by shear instabilities of outflowing very hot 
gas. The spectra of both of these phenomena are really produced by photoionization, but 
the photons are emitted from hot plasma instead of stellar photospheres. The contributions 
of metagalactic radiation upon the WIM at the largest distances from the plane may be 
important, especially for producing [O III]. While [O III] is surprisingly bright relative to Ha 
at large \z\ in NGC 891, the total power required for it is not large. The energy requirements 
of the WIM at modest \z\ ,say, ~ 1 kpc, are the most problematic constraint on the fraction 
of the WIM that is contributed by these alternative processes. 

The energy requirements of the extra heating can be determined from the details of the 
models. The detailed models show that the photoelectric heating from aT t = 50,000 K star 
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is 30% of the power required for the ionization. This number is approximately kT*/x(K)- 
This is expected, since fcT* is the mean energy of a black body spectrum above the ionization 
threshold x(H). For T* = 35,000 K, the extra heating power at \z\ = 1200 pc in the Perseus 
Arm (i.e., r_ 25 = 3) is 65% of the photoelectric heating, integrated over the volume, and 30% 
of the photoelectric heating for T* = 50,000 K. The power requirement of the extra heating 
is more modest in comparison to the total power needed, because we must then include the 
ionization energy in addition to the heating. At \z\ = 1200 pc the extra heating is 10% total 
power requirement of the WIM. The extra heating power, integrated over z from 500 pc to 
1200 pc, is ~ 5% of the total power requirement of the WIM. The power requirements for 
still larger \z\ are even smaller, since the intensity of the WIM decreases rapidly with height. 

While A3727 is the largest single coolant when extra heating is small, it is only about 
20% of the total cooling, with [N II] and [S III] not far behind it. The A3727 is, of course, 
very sensitive to extra heating because of its high excitation energy as compared to the 
other optical lines. Its importance relative to the other forbidden lines increases with T_ 2 5 
because its emissivity increases faster than that of the red lines. However, its overall share 
of the cooling remains about constant because of the onset of cooling by H° from collisional 
excitation of the Lyman lines. The collisional cooling exceeds A3727 at about T_ 2 5 = 1.5, but 
the precise boundary is set by the abundances and stellar temperatures. 

It is well known that [C II] 158 //m is a major coolant in the neutral ISM, as shown by 
models of photo-dissociation regions (PDRs) on the outer parts of molecular clouds (e.g., 
Kaufman et al. 1999). The COBE satellite showed (Wright et al. 1991) that 158 pim has 
a strength greater than the sum of all other infrared lines. The present models show that 
[C II] 158 /iin is not a major coolant in the WIM, although most of the gas phase C is 
singly ionized. The reason is that other coolants become dominant, while 158 fim is not 
boosted much by the relatively high temperature in the WIM as compared to PDRs. A 
simple calculation shows the effect. Let ji5&(T) be the emissivity (erg cm~ 3 s _1 ) of 158 fxm 
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and similarly for [N II] A6583. We find a mean T from setting Ixmss/Iua = J A6583<D / ju a (T) , 
where the j's involve the appropriate ionic abundances. The most favorable case for 158 /zm 
is using the largest (N/H), 75 ppm, and the lowest Ixm^/Inai 0.60, resulting in T = 7555 K. 
In this case, jim/ j&b&z = 0.33, which may seem significant. However, A6583 provides only a 
small part of the total cooling, according to models that predict A6583/Ho; = 0.6. For these 
models, the cooling of A6583 is ~12% of the total, so 158 /im is only 4% of the total. The 
Perseus Arm and NGC 891 models have lower (N/H) and (O/H) and, therefore, have hotter 
nebulae, making the 158 /xm even less important. 

The nebular temperature is raised by both extra heating and a shift of the contribution 
of hot stars to the ionization, so extra heating cannot be clearly identified until it dominates 
the photoionization heating. Figure 1 shows that this takes place for [O II]A3727 at T_ 2 5~ 
2. Below this value, the influences of stellar temperatures and abundances of coolants are 
comparable with those of extra heating. However, He + /H + is unaffected by extra heating. 
Very hot stars are needed if ^He approaches unity, especially for the low-density situation in 
the WIM. NGC 891 is an excellent example of this situation. 

It is significant that the combinations of models at low \z\ in the Milky Way cannot 
tolerate much extra heating, while extra heating is required at large \z\. 

The heating might arise from turbulence (Minter & Spangler 1997; Minter & Balzer 
1997). For models with uniform density, this extra heating can also be assumed to represent 
heating that is uniform per unit volume, such as magnetic reconnection (Raymond 1992; 
Gongalves et al. 1993; Birk, Lesch, & Neukirch 1998; Vishniac & Lazarian 1999). See Haffner 
et al. (1999) for discussion and references. 

Heating by PAHs is produced by a change of their state of ionization (Bakes & Tielens 
1994). Variations in the interstellar radiation field seen along the line of sight should affect 
the spectra of the near-infrared "Unidentified Infrared Bands" (UIBs) in the diffuse ISM and 
also the extra heating if PAHs are responsible for it. Observations of the spatial variations 
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of the UIBs, and correlations between these changes and the extra heating as shown by 
[N II] /Ha and [S II]/Ha, will be crucial in determining if PAHs are the main source of extra 
heating. Complete uniformity of the UIB spectra, suggesting that the distribution of stages of 
ionization of PAHs is uniform in various environments, would not support heating by PAHs. 

While there are many free parameters in the combinations of models, there are also 
strong constraints. The composition was taken to be the same at high and low \z\, and 
subject to observational values and the galactic abundance gradients for the local gas and 
the Perseus Arm. Up to 5 line ratios were fitted with 2 stellar temperatures and one relative 
fraction of ionizing photons. There are some combinations of observations that cannot be 
predicted by photoionization models. Our range of q considered (—4.80 < log(g) < —2.60) 
implies [N II] > 0.4Ha, with the minimum occurring for cool stars, T_ 2 5 = 0, and gas bounded 
models with large central voids. My models all have [O I] < 0.6H/3, with the maximum 
occurring for hot stars, the largest T_ 2 5 I considered (4.5), the lowest q, and a large central 
void. We never predict an [O III]A5007/H/3 < 0.15 combined with ^ Hc >0.9 (i.e., T* ~50,000 
K because of the generally low stage of ionization present, or small q). 

At the outer edge of the model, forbidden-line cooling (mainly from [S II]) cannot keep 
the temperature below ~ 10000 K if there is strong extra heating that is constant per H 
nucleus. Collisional ionization of H begins to produce a low but important fraction (> 
1%) of H + , along with corresponding electrons. In this case, large amounts of [S II] A6717 
and [O I]A6300 are produced (~ twice Ha) without any photoionization at all. Since such 
strong [O I] has not been observed, we conclude that collisional ionization and collisional 
excitation induced by the extra heating cannot produce most of the observed Ha. Clearly, 
a physical theory is needed to determine the level of H ionization at which various processes 
provide extra heating, and to estimate the conditions under which it becomes negligible in 
comparison to photoelectric heating. The present models arbitrarily end the extra heating 
when the fraction of H° is > 0.5, but values of 0.1 or 0.9 give essentially the same results. 
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However, the model cannot be continued to completely neutral H, which never occurs with 
constant extra heating because of collisional ionization of H. 

6. Summary 

As has been shown in past work, photoionization models can explain the observed 
spectra of the WIM. In this paper similar nebular models were applied to the WIM spectra 
from several rather different regions: the local spiral arm, the Perseus Arm at two values of 
\z\, and NGC 891 along a well observed strip perpendicular to the galactic plane. 

The combinations of nebular models assume that only two different stellar atmospheres 
provide the ionizing radiation. These stellar atmospheres (Kurucz 1991) are valid only in 
broad outline because they do not treat the stellar winds and non-LTE effects, and do 
not have the composition of the nebular gas. I regard the numerical values of the stellar 
temperature associated with each model with considerable caution because of these failings. 
However, WIM spectra basically count ionizing photons of various energies, so a distribution 
of stars with various temperatures can be simulated with two discrete values. Another 
question is whether the WIM is best characterized by combining stars of various temperatures 
and then computing the nebular spectrum, rather than combining nebular spectra from 
various individual stars. Other assumptions involve the geometry used in interpreting the 
model spectra, such as whether to use volume averages of integrals along various lines of sight 
through the model nebulae. Radiative transfer modeling of the fractal structure of the ISM 
is required before one can be confident about the propagation of the ionizing radiation from 
the source stars near the plane to large \z\. 

An important modification to the simplest picture of photoionizing the WIM is that 
extra heating beyond that supplied by photoionization is required (Reynolds et al. 1999). 
The diagnostics of this requirement are the large [N II] /Ha ratio that is observed at large 
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|^| in both the Perseus Arm and, especially, in NGC 891. Another feature is the observed 
constancy of [S II] /[N II]at ~0.6. There are several possible sources of this extra heating. It 
must be relatively unimportant for low \z\. An increase with density less rapid than n 2 helps 
fulfill this condition. 

Excepting NGC 891 at large \z\ (2 kpc in this paper), the models that are successful in 
fitting the WIM require both hot stars (T* ~ 45,000 - 50,000 K) and cool stars (T* ~ 35,000 
K), each providing comparable numbers of the ionizing photons. The cool stars are needed 
for the observed rather low values of [N II] Ha and [S II] /Ha. The hot stars are needed for 
[O III] and He+. 

NGC 891 at \z\ = 2 kpc is very interesting in that it shows very strong [O III] and He + , 
far above the plane. The combinations of models show that very hot stars (T* ~ 50,000 K) 
are needed, most plausibly at high \z\. If the exciting stars are found near the plane, their 
radiation must propagate through almost a vacuum up to the \z\ — 2 kpc region. 

The (N/H) required for the local gas is about what is observed from UV absorption 
studies, and that required for the Perseus Arm is about what is expected on the basis of 
local gas phase measurements plus galactic abundance gradients as observed in both stars 
and in H II regions. The fact that photoionization requires reasonable assumptions regarding 
the exciting spectra and compositions of the WIM makes it a very plausible candidate for 
explaining the principal physics of the WIM. 

With [N II] strengthened by extra heating, the [S II] follows if the (N/S) ratio is the 
same at high \z\ as near the plane. This constancy is expected if neither element is strongly 
bound to grains in dilute gas. The (S/N) required for the Perseus Arm is about 80% of solar. 

The abundances in both the Perseus Arm and the region sampled in NGC 891 seem to be 
lower in N, S, and O than the observed gas-phase abundances in the local arm. The accuracy 
of the composition predictions from these combinations of simple models is compromised by 



-26 - 



the model assumptions, but the result is so clear that it is almost surely real. 

This work has benefited from conversations with many colleagues, especially Ron 
Reynolds. Blair Savage, Bob Benjamin, Matt Haffner, and several colleagues at the workshop 
on the WIM held in Green Bank, WV, 23/24 September 1999 provided stimulating ideas and 
criticisms. I liked the suggestions from the anonymous referee. 
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Table 1 — Continued 
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a V>He = (He+/He)/(H+/H). 

b e is the one-cr accuracy with which the models are expected to reproduce the observations; 
see equation (1). 

c The prediction is with a combination of model H II regions with properties given in §[| in 
the text. 




Fig. 1— The strength of [0 II]A3727/H/3, plotted against the rate of extra heating in the 
models, for models with two different compositions and and various temperatures of exciting 
stars. The lines are all for the composition (C, N, O, Ne, S)/H = (85, 45, 250, 71, 8) ppM. 
The solid line shows T* = 50,000 K; the dotted, 40,000 K; the long dashed, 37,500 K; the long- 
short dashed, 35,000 K. The "error bars" indicate the maxima and minima for the models at 
various values of the excitation parameter, which determines the amount of H° at a particular 
fraction of the Stromgren radius. Those for cooler stars are displaced slightly to the right for 
clarity. The points show the same as the lines, except that the composition is (C, N, O, Ne, 
S)/H = (141, 75, 320,118, 13) ppM. The circles are for 50,000 K; the squares for 40,000 K; 
the triangles for 35,000 K. We see that extra heating strongly increases the strength of A3727. 
The composition, through changing the cooling, also enters the line ratio. 
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Fig. 2. — The same as Figure 1, except for the ratio [S III]A(9065+9531)/H<x We see that 
this line ratio is not a diagnostic of extra heating. 



